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Abstract 

It is shown that the free motion of massive particles moving in static space- 
times are given by the geodesics of an energy-dependent Riemannian metric 
on the spatial sections analogous to Jacobi’s metric in classical dynamics. In 
the massless limit Jacobi’s metric coincides with the energy independent Fer¬ 
mat or optical metric. For stationary metrics, it is known that the motion of 
massless particles is given by the geodesics of an energy independent Finslerian 
metric of Randers type. The motion of massive particles is governed by neither 
a Riemannian nor a Finslerian metric. The properies of the Jacobi metric for 
massive particles moving outside the horizon of a Schwarschild black hole are 
described. By constrast with the massless case, the Gaussian curvature of the 
equatorial sections is not always negative. 
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1 Introduction 


An elegant device for implementing the Principle of Least Action of Manpertnis was 
introdnced by Jacobi. One varies the action of a mechanical system 


Pidx'' = / piX^dt 


( 1 . 1 ) 


along an unparameterized path 7 in an n -dimensional conhguration space Q with 
coordinates x^ = 1 , 2 ,. ..,n and canonical momenta pi subject to the constraint 
that along the curve 7 the energy E is conserved. An equivalent formulation is to lift 
the curve 7 to the cotangent space T*Q and restrict variations to a level set of the 
Hamiltonian H{x,p) = E. 

In the simplest case the kinetic energy T = \mij{x)x^x^ , where the space depen¬ 
dent mass matrix mij{x)dx'^x^ endows the conhguration space Q with a Riemannian 
metric and the Lagrangian giving the equations of motion is 


L = -mij{x)xdx^ — V{x). 


( 1 . 2 ) 


Jacobi showed that the unparamaterised curves extremizing the constrained action 
are geodesics of the rescaled Jacobi metric 


jij{E, x)dx^dx^ = 2{E — V)mijdx'‘dx^ . 


(1.3) 


One recovers the parameterization of the motion as a function of physical time t by 
noting that length s with respect to the Jacobi metric is related to the original time 
parameter t by 


dt = 


ds 


(1.4) 


2{E-V) ■ 

This procedure opened up the way open to investigations of the motion of the 
original mechanical system using the methods developed by differential geometers to 
investigate geodesic motion. Of particular interest is the inhuence of the curvature of 
the Jacobi metric HIE]. An important application to gravity was the work by Ong 
[1] who studied the curvature of the the Jacobi metric for the Newtonian A^-body 
problem (see also 0 ). One has Q = with the flat metric 


mijdx^dx^ = madx^ 


2 

a 5 


( 1 . 5 ) 



and the potential energy is 


V = - T (1.6) 

Since in this case, the Jacobi metric is conformally flat the evalution of the curvature is 
straight forward. If = 2, the problem reduces to the Kepler’s problem of the relative 
motion and the relevant Jacobi metric is up to an unimportant over all constant factor 

a 


{— + —) (dr"^ + r‘^{d9‘^ + 9d(j)Y (1.7) 

m r 

By symmetry, one may restrict attention to the equatorial plane ^ = f which is a 
totally geodesic submanifold. One then has a 2 -dimensional axially symmetric metric 


{—+—){dr‘^ + r‘^d(p^) . ( 1 . 8 ) 

m r 

Ong pp showed that the sign of the Gaussian curvature of the metric has the opposite 
sign to that of the energy E. If E > 0 , which of course corresponds to unbound 
hyperbolic or parabolic orbits, he showed that the Jacobi metric fll.Sp is well dehned 
and complete for 0 < r < cxo and if E < 0 , which corresponds to bound elliptical 
orbits, it is well dehned for 0 < r < 

Ong also gave an isometric embedding of the Jacobi manifold into three dimen¬ 
sional Euclidean space with Cartesian coordinates x, y,z as a surface of revolution 
^ = /(\/x^ -|- y'^). If E = 0, then /" = 0 and the surface is the cone 2 ; = + y^ 

which has dehcit angle vr, or equivalently, semi-angle 30°. In the other two cases the 
surface approaches the cone near the origin. If E > 0 , then /" < 0 ; the surface has 
negative Gauss curvature and remains outside the cone . If E < 0, then /" > 0 the 
surface remains inside the cone and asymptotes the cylinder + y"^ = Ong 

also studied the three body problem using these techniques. 

It is obviously of interest, if only to extend one’s intuition by means of an easily 
visualised model, to see whether these ideas can be applied to General Relativity. At 
a formal level, one takes the conhguration space {Q, rriij} to be Wheeler’s superspace 
equipped with its DeWitt metric (cf. [1]). In the vacuum case, the potential V is 


Ry/gd? X 


(1.9) 


^Note that E in this Newtonian case does not contain a contribntion from the rest mass of the 


particle. 










and the Hamiltonian constraint implies that the energy vanishes i? = 0. One then 
obtains a picture of spacetime as a sheaf of geodesics in superspace. Having obtained 
the geodesic between two points in superspace one obtains the time duration between 
them using (ll.4p . thus solving the much discussed “problem of time”. 

Less ambitiously one may conhne attention to a mini-superspace truncation, and 
this has been done in attempts to investigate inflation [5] and the chaotic behaviour 
of Bianchi IX Mixmaster models (see e.g. 0 )- 

The focus of the present paper is different. It is the motion a test particle of rest 
mass m following a timelike geodesic in a stationary background. A limiting case 
would be a zero rest-mass particle. This latter case is well known in the static case 
to reduce to geodesic motion with respect to the optical or Fermat metric fij{x) = 
This was studied in [TJ [S] and a number of limitations on possible motions 
using the Gauss-Bonnet theorem. In particular in [H| the Gaussian curvature of 
the Schwarzschild optical metric restricted to the equstorial plane was shown to be 
everywhere negative and to approach a constant value near the horizon r = 2M. In 
[9] this behaviour was found to be universal for the near horizons of non-extreme 
static black holes. One purpose of the present paper is to extend this work to the 
case of massive particles. The extension of Fermat’s principle to cover stationary 
spacetimes entails replacing the Riemmanian metric by a Finsler metric of Randers 
type (see e.g. 


2 The Jacobi metric for static spacetimes 


If 


ds = —V'^dt^ -I- Qijdx^dx^ 


( 2 . 1 ) 


the action for a massive particle is 



( 2 . 2 ) 


where ^ The canonical momentum is 

at 

mx* 



(2.3) 







whence the Hamiltonian is 


H 


(2.4) 

(2.5) 


mV‘^ 

^ 1/2 _ g.-x^xi 


Setting 


Vi = diS, 

the Hamilton-Jacobi eqnation becomes 


+ V^gVdiSdjSj = E 


( 2 . 6 ) 


(2.7) 


or 


fWiSdjS = E^- , 

( 2 . 8 ) 

where fjk = Si and 



(2.9) 

is the optical or Fermat metric. Thus 



( 2 . 10 ) 


which is the Hamilton-Jacobi eqnation for geodesics of the Jacobi-metric jij given by 

jijdx^dx^ = — m‘^V‘^)V~'^gijdx'^dx^. ( 2 - 11 ) 


Note that the massles case, m = 0, the Jacobi metric coincides with the Fermat metric 
up to a factor of and as a consequence the geodesics, considered as unparameterized 
curves, do not depend upon the energy E. However in the massive case, m 7 ^ 0 the 
geodesics do depend upon E. 

In general, if the spacetime is asymptotically flat and the sources obey the energy 
conditions, then 0 < H <1. Therefore, if E"^ > m? , the Jacobi metric is positive 
dehnite and complete, even if horizons are present. If however E‘^ < there are 
bound orbits and the Jacobi metric changes signature at large distances. Generically 
there will be a level set of V on which E"^ — m?V‘^ vanishes and hence on which the 
Jacobi metric vanishes. From the point of view of the Jacobi metric this level set is a 
point-like conical singularity. Every geodesic must have a turning point on or inside 
this level set. 







3 The Schwarzschild Case 


In the case of the Schwarzschild solution the Jacobi metric is 

ds^ = (e" - + siu^^d^^^) . (3.1) 

The hrst large bracket in fl3.1l ) is the conformal factor and the second large bracket 
the optical metric. The later is dehned for 2M < r < oo and the horizon at r = 2M 
is inhnitely far way with respect to the radial optical radial distance or tortoise 
coordinate 


r* = 


' 2 M 


rf nr v * 

-^^ = r-2M + ln(—-1), r - 2M = 2MW{e^), (3.2) 


_ ZlVl 


'2M 


where W{x) is Lambert’s function dehned by Inx = W{x) + InkL(x). 

By spherical symmetry, in order to study geodesics, it is sufficient to consider the 
equatorial plane 6^ = | on which the restriction of the Jacobi metric is 


ds‘^ = ( + 


2 Mm? 


dr^ 


+ 


(X — (^X — — 


-defy 


(3.3) 


By axi-symmetry we have a conserved quantity often called Clairaut’s constant which 
corresponds physically to angular momentum. That is 


Now 


whence 


I = (E^ -m^ + 


2Mm? 


-m^ + 


2Mm? 


dr. 


d(j) 


(X - ^) ds 


+ 


ds' (l-2^)2 (X-Ml)Ms 


= constant. 




2 2 1 /dr.2 2M . 2 \ 

+—) (rv^(*) — )K+;i) . 


(3.4) 


(3,5) 


(3,6) 


which agrees with the standard result that 

.dr,n, . M, 




m' 




'dr' ' r ' ' 

where r is proper time along the particle’s worldline and 

2 ^4^ 


(3.7) 


I = mE 


dr 


(3.8) 




















as long as 


I- 2 M. 


dr = m- 


ds. 


(3,9) 


£2 _ „j2 + 2J£ml 

r 

Thus I is indeed angular momentum. In standard treatments u = ^ satisfies Binet’s 
equation 

(3.10) 


d^u F(u) 

+ u = 


d(j)‘^ ' "" h‘^u‘^ ’ 

where, for a massive particle orbiting a Schwarzschild black hole, we have 

F(u) 2 M 


(3.11) 


and where h = ^ is the conserved angular momentum per unit mass. If this were a 
classical central orbit problem we would say that we have a sum of an inverse fourth 
and inverse square law attraction. There is a first integral 

^ + C' = 2M(m - a)iu - I3){u - 7 ) (3.12) 

^d(p 


where C* is a constant related to the energy per unit mass 8 = — and the angular 
momentum per unit mass h by 


C = 


8^-1 


with 


a + /3 + 7 = 


2 M 


1 C 

P7 + 7« + ap = ^, 


(3.13) 


(3.14) 


3.1 Some Explicit Solutions 

The behaviour of the orbits depend on the two dimensionless quantities which are the 
specific energy 8 > 0 and In general, the solutions are given by elliptic integrals. 
However there is a one parameter family of explicit solutions of the form 

B 


u = A 


A = 


cosh"^ (cu^) 


6M 


(1±\1 


12MF 


1 / 12M2 

B = =f-^a/i 


2 M 


/l 2 


= ±-\ 1 


12 M2 


C = A\4MA-1). 


(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 




















These solutions arise because two of the the roots of the cubic eqns in fl3.12p coincide. 


1 


If oj is real and A + B ^ the solutions are symmetric about (/) = 0 , r = and 
end spiralling around a circular geodesic at r = ^. If = 12M^ then i? = 0 and 


A+B 


we have the inner most stable circular orbit at r = 6M for which the specihc energy 
£ = If h? = 16M^ the energy per unit mass £ = 1. Since A = —B = these 
orbits are in free fall from inhnity, starting from rest and spiral around a circular 
orbit at r = AM. All orbits starting from rest at inhnity (i.e. having £ = 1) with 
\h\ < AM fall through the horizon at r = 2M while all such orbits with \h\ > AM are 
scattererd back to inhnty. These latter orbits are relevant for the theory of the BSW 
effect fT3] . 

3.2 Bound States and Jacobi functions 

For a bound orbit we have three real positive roots taken to satisfy 


a>/9>n>7>0. 


(3.20) 


The hrst integral (13.121) leads to 



y/a — '-f du 


(3.21) 


y4(a - m)(/ 3 - u)(u - 7) 


with 



(3.22) 


Thus [16] 


u 



(3.23) 


where the modulus k of the elliptic functions is given by 



(3.24) 


and the quarter period K by 



2 


(3.25) 

















Using (I3.23P and properties of elliptic functions, (I3.2ip may expressed as 


1 1 I 1 ^ 1 I r 2f'^^-\\\ 

r =<—> +r" • ’ 


(3.26) 


where the constants Vp < Va are the radii at perihelion and aphelion respectively, since 
from fl3.20p 

(3.27) 

(3.28) 


r„ = — <r< — = r, 


and 


7 


^ ^ 2rpra 


/3 


rv + ra 


e = 


Tg - Tp 

Tg + Tg 


Perihelion is at : 0 = 

AK 

- n , n El 

w 

(3.29) 

Aphelion is at : 0 = 

AK ^ 1 

-(n + -), nez 

‘UJ 

(3.30) 


These expressions generalise the Newton-Kepler case for which the orbits are el¬ 
lipses with foci at the origin and are given by 


^ 7 cos2(|) + ^ sin2(|) = jil + e cos 0), (3.31) 

lip A f Q L Lj 

where L is the semi-latus rectum and e is the eccentricity Note that in both cases L 
is the harmonic mean of the perihelion and aphelion radii. 


3.3 Relation to Weierstrass Functions and Photon Orbits 

The term linear in u may be eliminated by setting 


U = V + c, 


(3.32) 


where 

2 c 1 


We then hnd that if 0 


c = 


1 

6M 


(i±\A 



1 — 6Mc 



5 


(3.33) 

0\/l — 6Mc, M = and C = C — c^(l — 2Mc), then 

(i) 2 + 2Mv^ + (7, (3.34) 

d(f) 


















which is the equation governing the orbit of a photon in a Schwarzschild solution of 
mass M. Thus the solution is [13] 


Mv = - + 2p{^), 
b 


(3.35) 


where p is the Weierstrass’s Elliptic function with parameter 93 = ^ — 
It follows that the massive particle orbits are given by 


1 — bMc 


u = 


p{(f)Vl — bMc). 


bM 3M 

A particular example is provided by the cardioidal photon orbit 

1 


V = 


2Mcos- 


2(h ’ 


which gives 


u=c+v=c+ 


1 — 6Mc 


2Mcos2(^dE|l£^; 


If we take the case for which 1 — bMc > 0 this gives 


(3.36) 
with (7 = 0, 

(3.37) 

(3.38) 


, 4M2 

bMu = 1 — a/ 1-——h 3- 


1 - 


4M2 


.(1-^)4^ 


(3.39) 


cos^ 


These orbits run from the past singularity at r = 0 out to a maximum radius 
given by 


bM 




1 + 2a/1 




(3.40) 


and return to the future singularity at r = 0. Note that r^ax > 2M. 

For a recent treatment of photon orbits in the Schwarzschild metric using Jacobi 
elliptic functions, the reader is directed to [T5] . 


3.4 Properties of the Jacobi meric 

Since any spherically symmetric is conformally flat we could adopt isotropic coordi¬ 
nates and avail ourselves of the results given in [1]. Alternatively, the calculation of 
the Gaussian curvature K with respect to the Jacobi metric on the equatorial planes 
is straight forward using equation (8) of [9]. However unless m = 0, this leads to 
rather complicated and un-illuminating expresions. A qualititive analysis, based on 
















the behaviour of circular geodesics, i.e. those with r = constant seems preferable. 
We restrict attenton to an equatorial plane. 

li > m?, K tends zero at inhnity and in all three cases it tends to — 
as r tends to the horizon at r = 2M, which is at inhnite Jacobi radial distance, the 
curvature tends to a negative constant. If E'^ < m?, the Jacobi manifold has an outer 
boundary at which the metric vanishes. This happens when 

2M. 2Mm2 


E^ = m\l 


r = 




h;2 


> 2M. 


(3.41) 


This outer boundary should be thought of as a point since the Jacobi circumference 


2M, 


27r{E‘^ -m^{l -)) 


(3.42) 


— ?M. 

r 

vanishes there. In the vicinity of the boundary K is positive. 

Circular Jacobi geodesics are possible. These correspond to extrema of the Jacobi 
circumference and are located at values of r for which 
3M, 




2M.2_q E^_{r-2Mf 

r ’ m2 r(r — 2M) 


(3.43) 


Circular geodsics exist with real energies of all values of r > 3M. For every value of 
^ > 1 there is a unique circular geodesic with radius r between 3M and AM. For 
every value of ^ between | and 1 there are two circular null geodesics, the inner, 
which is unstable, has its radius between AM and 6M, and the outer whose radius is 
greater than 6M. Three interesting cases arise 

• m^ = 0 , r = 3M. These are circular null geodsics which are cirular geodsics 
of the optical metric. 

• E^ = rn? , r = 4M, these are have the smalles raduius among all circular 
bound geodesics. 


^2 = |m2 


r = 6M. These are the most deeply bound circular timelike 

geodesics. 

These results are consistent with the standard approach to circular timelike geodesics 
which is to require the simultaneous vanishing of the right hand side of (13.71) and its 
derivative. Eliminating E'^ and solving for u gives 

1 


6M 


12M2\ 


h2 


-}■ 


ru = 


1± A 1- 


(3.44) 
















Elimininating E and solving for h gives 


16(r — 3M) 

-■ ( 3 -«) 

For the unbound circular geodesics with E'^ > m? , varies from 0 to | as r 

varies from 3M to AM. For the bound circular geodsics with w? < E^ < |m^ one 
hnds that varies from | to 1 as r varies from AM to 6M where it achieves its 

maximum value and therafter as r varies from 6M to inhnity it decreases monotoni- 
cally to zero. 


3.5 Gauss Curvature and Isometric Embedding 

We can use these results to say something about the Gauss curvature K. 

We begin by recalling that if the induced metric on the surface of revolution 
^ = fip) > P = \/x2 +|/2^is 

A^dr"^ + = (1 + {f'Y)dp^ + p^dcj)^ (3.46) 


then we have 


1 dr 

p = G(r), (/')-! = (-—)2-l). (3.47) 

In our case, the embedding will extend from inhnity towards the horizon at r = 2M 
as long as the r.h.s of fl3.47p remains positive. The Jacobi metric is conformal to the 
spatial metric of the Schwarzschild solution and its equatorial plane is well known 
to be isometrically embeddable as all the way down to the horizon as the Flamm 
paraboloid [Ml HH] 


z = 


{p-2Mf 

8M 


(3.48) 


However this is in general not possible for the Jacobi metric. 

For example if m = 0 one is limited us to the region r > |M. In that case 


K 


2M 3M 


(3.49) 


which is everywhere negative and near the horizon the optical metric is asymptotic 
to one of constant negative curvature equal to 

This is analogous to the well known fact that the metric with A = 1 and C = 
ae~^ , a > 0 has constant negative curvature — ^ and may be embedded into as the 













surface of revolution whose meridional curve is a tractrix. However this “ Beltrami’s 
trumpet” for which 


/' = --v/a^ — , ±/= constant + a cosh ^(-) — y/a? — ^ (3.50) 

P P 

is incomplete and only occupies the region for which p < a, i.e. r > 0. 

By fl3.3p the Jacobi metric is conformal to the optical metric with conformal factor 

2Mrn^ 


E^-m^ + 


(3.51) 


It also tends to constant negative curvature near the horizon. Thus in general one does 
not expect to be able to capture the near horizon geometry of the Jacobi metric by 
an isometric embedding into as a surface of revolution. This represents a practical 
obstruction to constructing black hole analogues using such materials as graphene 


Consider now the case > |m^. Near the outer boundary the Gauss 

curvature is positive and it is positive at the outer circular geodsic which is a local 
maximum of the Jacobi circumference. By the time we get to the inner, unstable orbit, 
which is a local minimum of the Jacobi circumference, the Gauss curvature is negative 
and it is negative near the horizon at r = 2M near which the Jacobi circumference 
diverges. If E"^ < |m^ there are no circular geodesics and the curvature is positive 
near the boundary and negative near the horizon. Thus the Gauss curvature of the 
Jacobi-metric restricted to the equatorial plane is not everywhere negative as is the 
case for the Fermat metric. 


4 The Jacobi metric for stationary spacetimes 


We cast the spacetime metric in Zermelo form HU 

112 


g^ydx^dx'’' = 


1 - hijW^W^ 


dE + hij{dx^ - W^){dx^ - W^) 


(4.1) 


where I shall call hj, the Zermelo metric, IT* the wind, and 


= -gu , 


(4.2) 


where = ^, is the timelike Killing vector held. Note that if the wind vanishes, 

the Zermelo metric coincides with the optical or Fermat metric. 








The Lagrangian L for a point particle of mass m undergoing geodesic motion in 
a spacetime with metric fl4.ll) is 

mV 


L = 


- hijWW 


hij{V - W^){x^ - W^) (4.3) 


where V = The canonical momenta pi are therefore given by 


Pi = 


mV 


hij{x^ — W^) 


- h^jW^W^ - hij{xi - W^)(x^ - W^) 
The Hamiltonian, H = piV — L, is given by 


m^ 


H{m, Pi, x*) = h^^piPj + —1 - hijWW^ + PiW ^, 


(4.4) 


(4.5) 


where h^^hkj = 5). Note that in the massles case, m = 0 the Hamiltonian H becomes 


H{0,pi,x'-) = ^ h^^PiPj +PiWi 


(4.6) 


which coincides with equation (14) of [10] thus recovering the result that the pro¬ 
jection of null geodesics of a stationary spacetime onto the space of orbits of the 
timelike Killing vector solve the Zermelo problem of minimizing the time of travel in 
the presence of the wind and with respect to the Zermelo metric hij. 

A quick way of obtaining (14.51) is to work on the co-tangent bundle of spacetime. 
The so-called super-Hamiltonian "H is subject to the constaint 

= g^'^Pi^Pu = -m^ • (4.7) 


If one solves 04.71) for po = —H one obtains 04.5p . 

In general, on the level set H = E of the Hamiltonian we have 


h'^PiPj + 


(1 - hijW^W^) 


= {E-piWy. 


(4.8) 


If the mass m is non-zero, it is not possible to cast this in he form of an expression 
which is a homogeneous degree two in momenta pi equated to a constant. If it were 
so, then and a Legendre transform would result in a Lagrangian which is of degree two 
in velocities and hence we would be dealing with a Finsler structure, possibly 
Riemannian, as in the case of a static metric. Thus we are faced with a geometric 
structure more general than a Riemannian or even a Finsler metric. 
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